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IS,  ABSTRACT 


The  Maximum  Principle  was  applied  to  determine  the  types  of  temperature  variations 
that  minimize  and  maximize  the  breadth  of  the  molecular  weight  distribution  (MWD)  for 
chain  addition  polymerizations  in  batch  reactors.  It  was  found  that  the  variations 
which  minimize  the  breadth  of  the  MWD  keep  the  instantaneous  number  average  chain 
length  constant.  -The  variations  which  maximize  the  breadth  of  the  MWD  are  step  changes 
in  temperature  resulting  in  bimodal  distributions.  Numerical  and  experimental  examples 
of  such  variations  are  presented.  MWDs  with  minimum  and  maximum  breadths  are  compared 
to  those  that  might  be  formed  by  temperature  variations  in  real  reactors.  Under  most 
conditions,  temperature  variations  appear  to  have  a  much  greater  effect  on  MWD  than 
residence  time  distributions  and  mi cromi x i ng. 
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ABSTRACT 


The  Maximum  Principle  was  applied  to  determine  the  types  of  temperature 
variations  that  minimize  and  maximize  the  breadth  of  the  molecular  weight 
distribution  (MV/D’)  for  chain  addition  polymerizations  in  batch  reactors.  It 
was  found  that  the  variations  which  minimize  the  breadth  of  the  MWD  keep  the 
instantaneous  number  average  chain  length  constant.  The  variations  which 
maximize  the  breadth  of  the  MV/D  are  step  changes  in  temperature  resulting  in 
bimodal  distributions.  Numerical  and  experimental  examples  of  such  variations 
are  presented.  MWD's  wi th  minimum  and  maximum  breadths  are  compared  to  those 
that  might  be  formed  by  temperature  variations  in  real  reactors.  Under  most 
conditions,  temperature  variations  appear  to  have  a  much  greater  effect  on  MV/D 
than  residence  time  distributions  and  micromixing. 

4 

INTRODUCTION 

Chain  addition  polymers  are  formed  by  the  addition  of  monomer  units,  one 
at  a  time,  to  growing  chains.  These  chains  have  life  times  which  are  much 
shorter  than  the  half  life  of  the  monomer.  Polymers  produced  by  free  radical 
mechanisms,  such  as  polystyrene,  low  density  polyethylene,  poly-methyl  methacrylate, 
etc.,  are  examples  of  chain  addition  polymers. 

Denbigh  (1]  first  reported  that  polymers  produced  in  commercial  reactors 
would  not  generally  have  the  same  molecular  weight  distributions  (MWD),  and 
thus  physical  properties,  as  polymers  produced  in  laboratory  batch  reactors. 

In  particular  he  investigated  the  effect  of  residence  time  distributions  (RTD) 
on  MWD.  Tadmor  and  Biesenberger  [2]  quantified  the  RTD  effect  and  investigated 
the  effect  of  micromixing  as  well.  They  found  that  the  effects  of  RTD  and 
micromixing  on  MWD  were  relatively  small  for  chain  addition  polymers.  These 
results  were  in  agreement  with  the  qualitative  predictions  of  Denbigh,  who 
reasoned  that  since  the  life  times  of  growing  chains  were  short  compared  to 


the  average  residence  time  in  a  real  ■'eactor,  RTD  should  have  little  effect 
on  MWD.  However,  the  important  e^.e-.ts  of  temperature  variations  have  not 
been  studied  in  such  detail.  Cr.<  !  addition  polymerizations  are  highly  exothermic 
AH  =  -15  to  -25  kcal/gmol.  This  'v.upled  with  the  high  viscosity  and  low  thermal 
conductivity  of  polymer  system;,  akes  heat  transfei  and  thus  temperature  control 
difficult.  In  fact  non-isothcr  vl  operation  of  bulk  polymerization  reactors  is 
probably  the  rule  rather  than  . e  exception. 

The  problem  of  determini  \  the  effect  of  temperature  variations  on  the 
MWD  can  be  attacked  in  two  w  /■  .  The  first  is  to  model  specific  reaction  systems. 
To  some  extent  *h!s  has  beei  d'-ne  for  tubular  reactors  [3],  [*♦] ,  [5]  •  The  second 
way  is  to  try  .o  bound  the  • <ects  of  temperature,  i.e.,  to  determine  what 
temperature  variations  mir  m'ze  or  maximize  the  breadth  of  tne  MWD.  The  second 
approach  has  been  chosen  for  this  study.  This  approach  has  the  advantage  that 
optimal  control  theory  can  ba  applied,  and  leads  to  generalized  results. 

Hoffman,  et  al.  [f>J  and  others  [7.8,9],  from  a  kinetic  understanding  oF 

the  polymerization  process,  proposed  temperature  and  monomer  or  initiator  addition 

policies  which  minimize  the  breadth  of  the  MWD.  These  policies  keep  the  insan- 

taneous  average  chain  length,  X"  .,  constant  throughout  the  reaction.  Howevt”, 

n » i 

they  are  not  always  physically  realizable  because  of  constraints  on  the  reditant 
concentrations  and  on  the  temperature.  Ray  and  his  co-workers  [9]  made  nume  ical 
studies  for  a  styrene  polymerization  model  using  optimal  control  theory.  The 
best  sub-optimal  policies  were  found  to  be  very  close  to  the  policies  mentioned 
above. 

In  this  paper  the  general  nature  of  the  temperature  variations  that  minimize 
and  maximize  the  breadth  of  the  MWD  will  be  deduced  from  the  Maximum  Principle. 

The  polymerization  model  to  be  considered  takes  the  gel  effect  into  account  emperi 
cally.  Theoretical  results  are  compared  with  experiments.  Also,  temperature 
policies  which  might  exist  In  real  reactors  are  discussed  comparatively . 


I 
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KINETIC  MECHANISM 

A  typical  free  radical  mechanism  [10]  with  chemical  initiation  was  chosen 

t  h 

to  describe  the  kinetics  of  the  system.  This  mechanism  consists  of  n  order 
initiator  decay,  first  order  proportion  with  respect  to  monomer,  and  second 
order  termination  by  combination  and  disproportionation.  The  absence  of  transfer 
reactions,  and  the  principle  of  equal  reactivity  of  growing  chains  were  assumed. 
The  mechanism  can  be  represented  by  the  following  reactions: 


It  is  well  known  that  many  chain  addition  polymerizations  are  auto-catalytic 
[10].  This  phenomenon  has  been  termed  the  gel  effect,  and  abundant  experimental 
evidence  indicates  that  It  is  due  to  a  diffusion  controlled  termination.  Theore¬ 
tical  [11]  and  emperica!  [12]  relationships  between  kt  and  solution  viscosity  have 
been  proposed  t  account  for  the  gel  efrect.  Since  the  viscosity  of  polymer 
solutions  is  primarily  dependent  on  polymer  concentration,  as  a  first  approximation, 
the  overall  termination  constant  was  assumed  to  be  an  empericalj separable  function 
of  monomer  conversion,  m,  and  reciprocal  temperature,  y  =  1/R^T,  to  account  for 
the  gel  effect. 

kt  =  ktc  +  ktd  =  kto(y)-  g2^  '  At  exp  (*Ety^‘  g2^  ^ 

Assuming  constant  density  and  that  primary  radicals  do  not  terminate,  and 
applying  the  quasi  steady  state  approximation  [  1 3 ]  ,  the  following  equations  for 
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the  conversion  of  monomer  and  initiator  and  for  the  first  three  moments  of 


the  MWD  in  a  batch  reactor  can  be  derived. 

5T  =  ft4  =  k|  ly)  ■  lr  0‘c)  0-»)/S  +  2ak2(y).  h(1-c)  (7) 

=  k2 (y) •  h(l-c)  (8) 

=  (2-n)  ak2(y).  h ( 1  -c)  (9) 

Jjf*-  2(l+n)  k,(y).h40-c)  (l-m)/g  +  (2+n)k3(y). 

(1-m)  +  2riak2(y)*  h(l-c)  (10) 


Where  h(l-c)  is  any  arbitrary  function  of  initiator  conversion. 
Details  of  the  derivation  of  eqs.  7“ 1 0  are  given  in  Appendix  A. 
The  number  and  weight  average  chain  lengths  are  then  given  by: 


m/Co  =  m/((2-q)ac] 

(ID 

(C2  +  m)/m 

(12) 

The  instantaneous  number  average  chain  length,  i.e.,  the  average  lergth 
of  dead  chains  being  formed  at  any  instant,  is  given  by: 

7  dm  J _  .  kl(y)  .  (l-m)  +  2 

n,  i  dCo  =  (2-f|)a  k2  (y)  hi  (l-c)g  (2-n)  03) 

Eq.  13  shows  that  the  gel  effect,  decreasing  g,  will  result  in  an  increase  in 

the  average  chain  length.  For  typical  chain  addition  polymerizations  E2  i  30 

kcal/mol .  while  Ej  -  20  koal/mol.  Thus  kj/k2  has  a  negative  activation  energy, 

and  I  .  varies  inversely  with  temperature.  Finally,  two  modes  of  polymerization 
n,  i 

may  te  denoted  from  eq.  13  depending  on  the  relative  rates  of  conversion  of 
monomer  and  initiator.  In  'conventional'  polymerizat ion  the  rate  of  initiator 
conversion  is  slower  than  the  rate  of  monomer  conversion.  X  .  decreases  with 
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tinie.  In  many  'conventional'  polymerization  systems  the  rate  of  initiator 

conversion  is  so  low  that,  with  negligible  error,,  the  initiator  concentration 

may  be  assumed  to  be  constant  throughout  the  run.  In  'dead  end'  polymerization 

[ 1  ]  the  rate  of  initiator  conversion  is  faster  than  the  rate  of  monomer  con- 

>  .rsion.  X  .  increases  with  time, 
i*.  i 

MINIMUM  BREADTH 

The  objective  of  the  problem  is  to  determine  what  types  of  temperature 
variations  result  in  the  narrowest  MWD's.  The  dispersion  index: 

0  -  X/X  =  1  +  (o/X  )2  (14) 

n  w  n  n 

will  be  used  as  a  measure  of  the  breadth  of  the  MWD.  To  facilitate  comparing 
results,  Dn  will  be  minimized  at  some  predetermined  monomer  conversion,  m*, 
and  number  average  chain  length,  X  *.  The  temperature,  and  thus  the  reciprocal 
temperature,  y,  will  be  assumed  to  be  physically  constrained 

T*  <  T  *  T*;  y*  *  y  *  y*  (15) 

These  temperature  bounds  might  be  the  temperature  of  the  cooling  water  in  the 
reactor  jacket  and  the  ceiling  temperature  of  the  polymer.  Using  eqs.  11  and  12, 
the  problem  can  be  formally  stated  as  follows  for  the  case  of  constant  n:  Find 
the  temperature  policy,  y(t),  subject  to  the  constraints  y*  f  y  $  y^,  that  minimizes 
$2  (0)  at  a  given  m*(f>)  and  c*{0). 

For  constant  n  the  state  of  the  system  is  described  by  eqs  7,  8,  and  10. 

The  Maximum  Principle  states  that  the  optimal  policy,  y(t),  which  meets  the 
problem  objective,  maximizes  the  Hamiltonian  and  makes  it  equal  zero.  Since 
we  are  considering  only  one  control  variable,  y,  these  criteria  can  be  formalized 
to  yield: 

H  *=  0 

*  °;  <  o  ;  y*$yEy,$y* 

ay2 


when  y  is  unconstrained: 

5H 

3y 


(16) 


when  y  is  at  a  constraint 


3H  >0 

! 

*< 

n 

3y 

3H  <0 

> 

y  -  y* 

ay 

The  Hamiltonian  and  adjoint  variables  are  derived  in  Appendix  B.  Eq.  16 
can  then  be  applied  to  obtain  the  criteria  for  the  optimal  policy.  Details 
of  this  procedure  are  also  in  Appendix  B.  When  n  and  f  are  constant,  the  re¬ 
sulting  criteria  for  the  optimal  policy  are: 


kj(l-m)  + 

2  <  U,-2-2n)  +  2  .  y  «  y> 

X  • 
n,i 

c  ■  ■  ■  ■  — - 

(2-n)ak2h^g 

(2-n)  (2+n)  (2-n)  (2-n) 

(17a) 

7  . 

=  k,  (1-m) 

+  2  const  ;  y*  $  y  =  y$S  y* 

(17b) 

n,  i 

*  i 

(2-n)  ak2h‘g 

(2-n) 

r  , 

n,  1 

„  (1-m) 

♦  2  >  <V2-2">  *  2  ; 

y  -  y* 

(2-n)ak2h^g 

(2-n)  (2+n)  (2-n)  (2-n) 

(17c) 

The  stationary  policy  (17b),  i.e.,  the  policy  that  makes  3H/3y  “  0,  is 

identical  to  the  policy,  postulated  from  a  purely  physical  argument  16,7,8,9], 

which  keeps  X  .  constant, 
n  9 1 

Physically,  the  optimal  policy  appears  to  minimize  Dn  by  minimizing  the 
change  in  the  MWD  caused  by  changing  monomer  and  initiator  concentrations  during 
polymerization.  Referring  to  eq.  13,  if  j  can  be  considered  to  be  the  product 
of  a  temperature  dependent  term,  kj/k2>  and  a  concentration  dependent  term, 

(l-m)  .  The  concentration  dependent  term  is  proportional  to  the  isothermal 

(2-n)a  h*g 

>T  ,  where  the  isothermal  X  .  is  calculated  at  a  temperature  that  results 
n,  i  n» 1 

In  the  same  m*  and  Xn*  as  the  optimal  policy.  Three  possible  physical  cases 

exist:  I)  isothermal  X  .  decreases  wi th  time,  'conventional  polymerization'; 

n ,  i 
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II)  isothermal  Xn  .  increases  with  time,  'dead  end  polymerization'  or  very 

strong  gel  effect;  and  III)  isothermal  X  .  decreases  then  increases  with 

n  ,  t 

time,  moderate  gel  effect.  For  each  case  three  sub-cases  exist  depending 
on  the  initial  temperature  of  the  system.  The  initial  temperature  is  in  turn 
a  function  of  the  target  m*  and  Xn*,  the  initial  concentration  of  reactants, 
and  the  initiator  employed.  The  nine  cases  are  summarized  below: 

I.  Isothermal  X  .  decreases  with  time. 

n,  i 

a.  y(0)  =  yA,  Inequality  17a  insures  that  the  temperature 
remains  equal  to  yv.  throughout. 

b.  y(0)  =  ys;  Temperature  decreases  to  keep  X^  .  constant 
and  may  reach  yA,  upon  which  the  policy  switches  from  y&  to  yA.  Inequality 

17a  insures  no  further  temperature  change  occurs. 

/  \  *  * 

c.  y (0)  -  y  ;  Temperature  remains  y  until  the  Inequality 

17c  comes  to  equality.  When  and  if  this  occurs,  the  policy  switches  to  ys, 
a  decreasing  temperature  policy.  If  y$  decreases  to  y^.,  the  policy  switches 
to  Yj..  Inequality  17a  insures  that  no  further  temperature  change  occurs. 

II.  Isothermal  X  .  increases  wi^h  time. 

n ,  i 

5.  y(0)  «  y*;  Temperature  remains  yA  until  Inequality  17a 

comes  to  equality.  When  and  if  this  occurs,  the  policy  switches  to  ys, 

an  Increasing  temperature  policy  to  keep  X  ,  constant.  If  y  reaches  y*, 

n  1 1  s 

the  policy  switches  to  y*.  Inequality  17c  insures  that  no  further  temperature 
change  occurs. 

b.  y (0)  «=  y  ;  Temperature  increases  to  keep  X  ,  constant. 

s  n ,  i 

If  ys  reaches  y*,  the  policy  switches  to  y*.  Inequality  17c  Insures  that 
no  further  temperature  change  occurs. 

c.  y(0)  =  y*;  Inequality  17c  insures  that  the  tr nperature 
remains  y*  throughout. 

III.  Isothermal  X  .  decreases  then  increases  with  time. 

n ,  i 

a.  y(0)  *  y*  Same  as  I  la 

b.  y(0)  =  y  On  the  decreasing  portion  of  isothermal  7  .  (t) 

5  n  .  i 


I 


I 

*  .  ! 

.  '  3 

i 

i 

the  temperature  decreases  to  keep  Xn  .  constant  and  may  reach  y^,. ,  where  the 

policy  switches  to  y.,_.  ,0n  the  increasing  portion  of  isothermal  X*^  .  (t)  the 

temperature  increases  to  keep  X  .  constant,  and  may  reach  y*  where  the  policy 

n,  i 

switches; to  y*.  Inequality  17c  insures  that  no  further  changes  are  possible. 

I  < 

c.  y(0)  =  y*.  On  the  decreasing  portion  of  isothermal 

Xn  ;  (t)  inequality  17c  may  come  to 'equality  and  the  policy  switches  to  ys> 

a  decreasing  temperature  policy  and  may  reach  y..;.  If  the  switch  from  to  ys 

occurs,  the  temperature  increases  on  the  increasing  portion  of  isothermal  X  (0 

n ,  i 

and  may  reach  y*  again.  If  the  switch  from  y*  to  yg  does  not  occu*,  the  policy 
is  y*  throughout.  i 

i 

Maximum  8readth 

The  objective  of  this  problem  is  to  determine  the  types  of  temperature 
variations  that 'result  in  the  broadest  MWD's.  The  restrictions  imposed  in 
the  minimum  Dn  case  will  again  apply.  The  problem  may  be  formally  stated 
as  follows:  Find  the  temperature  policy,  y(t),  subject  to  the  constraints 
y*J  -$y  that  maximises  (0) . 

for .constant  n  the  state  of  the  system  is  described  by  eqs.  7,  8  and  10. 

The  Hamiltonian  and  adjoint  Variables  are  derived  and  eqs.  16  are  applied  to 

I 

obtain  the  criteria'for  the  optimal  policy.  Details  of  this  procedure  can 
be, found  in  Appendix  C.  When  f  is  constant,  the  resulting  critera  for  the 
optimal  pol icy  are: 


X  ?  -  _ 

'•  (l-m) 

1 

2 

t- 

<  .  (X,+2+2n)  + 

? 

;  y=y* 

n,  i 

(:>.*n;ak2h^g 

(2-n) 

(2+  n)  (2-n) 

(2-n) 

(18a) 

Xn5«_ 

k,  O-m) 

+  2 

>  .  (X j+2+2r,  )  + 

2 

.  ;  y=y* 

n ,  j 

(?*r,)ak2h^g 

(2-n  ) 

(2+n  )  (2-n) 

(2-n) 

(18b) 

The  stationary  policy,  y*  <y^,  cannot  *'or.i,  part  of  the  optimal  path.  The 
nolicy  that  maximizes  Dn  exists  entirely  on  the  bounding  temperatures,  y*  ai 
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,  and  therefore,  the  optimal  temperature  policy  must  consist  of  y:':  or  y... 
or  step  changes  between  the  two. 

The  number  and  direction  of  changes  betv/ecn  y*  arid  yA  is  determined  by 
the  sign  of  Aj,  the  nature  of  the  change  in  the  isothermal  X^  .  with  time, 
and  the  initial  temperature  of  the  system. 

If  Aj  is  positive,  18  states  that  the  temperature  must  be  y#.  for  any 
initial  conditions  and  target  m*  and  X^*.  This  is  not  the  usual  case  since 
a  policy  of  y  =  y^.  will  not  generally  produce  a  polymer  with  the  target  conditions, 
m*  and  Xn*.  Therefore,  Aj  will  generally  be  negative.  When  Aj  is  negative, 
six  cases  can  be  identified: 

I.  Isothermal  X  .  decreases  with  time: 

n,  i 

a.  y(0)  =  y,. .  Inequality  1 8b  may  come  to  equality  and  the  policy 
switches  to  y*.  Further  switches  are  precluded  by  1 8a . 

b.  y(0)  =  y*.  Policy  is  y*  throughout. 

II.  Isothermal  X  ,  increases  with  time. 

n.i 

a.  y(0)  =  y^.  Policy  is  y.v  throughout. 

b.  y(0)  =  y*.  Inequality  1 8a  may  come  to  equality  and  the  policy 
switches  to  yv. .  Further  switches  are  precluded  by  18b. 

III.  Isothermal  X  .  decreases  then  increases  with  time. 

n,  i 

a.  y(0)  =  y*.  Inequality  1 8b  may  come  to  equality  on  decreasing 

segment  of  X^  .  vs.  t  and  policy  switches  to  y*.  If  this  occurs,  inequality 

18a  may  come  to  equality  on  the  increasing  segment  of  X  .  vs.  t,  and  the  policy 

n,  i 

switches  to  yA.  Further  switches  are  precluded  by  18b. 

b.  y(0)  =  y*.  Inequality  1 8a  may  come  to  equality  on  the  in¬ 
creasing  segment  of  X  .  vs.  t  and  the  policy  switches  to  v^.  Further  switches 
are  precluded  by  1 8b . 

Since  an  isothermal  policy  of  y*  or  yft  will  not  generally  reach  a  pre¬ 
determined  m*  and  X  *,  most  policies  that  maximize  D  will  consist  of  one  or 
n  K  n 

at  most  two  step  changes  in  temperature,  Cases  la,  Mb  and  jllaSb.  'Jnlike 


the  minimum  On  cose,  the  physical  stategy  behind  the  optimal  criteria  is  not 
evident  from  18.  However,  this  stategy  will  become  evident  when  numerical 
examples  are  considered, 

ISOTHERMAL  X  , 

- n.i 

In  the  preceedlng  sections  the  general  nature  of  the  optimal  policies 
that  minimize  and  maximize  the  breadth  of  the  MWD  was  found  to  depend  on 
the  variation  of  the  isothermal  Xn  .  with  time.  This  variation  may  be  quanti¬ 
fied  by  studying  dX  ,/dt  at  any  fixed  temperature  level. 

n  y  1 

Differentiating  eq.  13  with  respect  to  time  and  substituting  eq.  7  and  8 

into  the  resui t.yields  dX  ./dt  at  any  fixed  temperature  level  as  a  function 

n ,  i 

of  m  and  c.  Dividing  eq.  8  into  eq.  7  after  making  the  long  chain  approximation 
[10],  neglecting  the  second  term  on  the  right  h-.id  side  of  eq.  7,  and  integrating, 
yields  c  as  a  function  of  m  in  an  Isothermal  system.  With  this  result,  dX^  ,/dt 
can  be  rearranged  to  give  the  following  criteria: 


(k2/2k,) 

.  k2  .W(m) 
2kl 


g  +  (l-m)g'  <  0  ;  Isothermal  X  .  decreases 

-  (19) 

g2  >  0  :  Isothermal  X  ,  increases 

s  n,i 


where  W(m)  =  /  g  lit/)  dm*/ ( 1 -m)  and  g'  =  dg/dm.  Eq.  19  was  used  to  prepare  maps 

*'o 

of  kj/kj  var'ous  temperature  levels  versus  conversion  that  show  regions 

where  isothermal  X"  .  increases  and  decreases.  Maps  for  styrene,  methyl  metha- 
n,  i 

crylate  and  for  polymerization  without  the  gel  effect  are  shown  on  Figure  1. 


CONTROLLABILITY 

In  the  preceedi ng  sections  it  was  tacitly  assumed  that  it  was  possible 
to  reach  the  target  conditions,  m*  and  X^*,  by  some  constrained  temperature 
policy,  Ta<  T$  T*.  It  may  not  always  be  possible  to  do  this,  however,  because 
of  insufficient  reactant,  concentrations  or  temperature  constraints  that  are 
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too  narrow.  Such  systems  are  not  controllable. 
The  limits  of  control labi 1 iby  are  given  by: 


n  v 
T=T* 


X  * 
n 


m=rii'" 

T=Tj. 


a 


m*  1_ 

TF^'y.- 


(20) 


(21) 


NUMERICAL  EXAMPLES 

Two  vinyl  polymerization  systems,  poly-styrenc  and  poly-methyl  methacrylate 
initiated  by  aso-bis-iso-butrylni tri le  (A1BN),  were  chosen,  and  numerical  examples 
of  optimal  temperature  policies  were  calculated. 

The  initiator  decay  is  first  order,  h(l-c)  =  1-c. 

The  kinetic  data  that  were  used  are  summerized  in  Table  1. 

Kinetic  constants  for  styrene  were  taken  from  Hamielec,  et  al.  [15].  The  gel 
effect  function,  g,  was  obtained  from  the  data  of  Tobolsky,  et  al.  [16]  and 
Nishimura  [17]-  Bevington,  et  al.  [18]  showed  that  termination  for  polystyrene 
occurs  exclusively  by  combination,  n=l.  Kinetic  constants  for  methyl  methacrylate 
were  chosen  by  comparing  the  individual  rate  constants  determined  by  various 
workers  [18,19,20,21,22]  to  the  apparent  rate  constant,  kt/kp  ,  obtained  by 
Bevington,  et  al.  [23]  and  by  Baysal  and  Tobolsky  [ 2^ ] .  On  this  basis  the  in¬ 
dividual  rate  constants  used  by  King  and  Skaates  [25]  were  chosen.  The  gel  effect 
function  was  obtained  from  the  data  of  Hayden  and  Melville  [19].  The  data  of 
Bevington,  et  al.  [18]  and  of  Obrian  and  Gornick  [26]  for  methyl  methacrylate 
show  that  n  is  a  function  of  temperature.  However,  the  primary  mode  at  temperatures 
of  interest  appears  to  be  disproportionation,  n=0,  and,  therefore,  termination 
was  assumed  to  occur  exclusively  by  this  mechanism. 

Optimal  policies  were  calculated  on  a  POP- 1 0  digital  computer.  The  cal¬ 
culation  procedure  used  was  as  follows: 
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1 .  Guess  X | (0) 

2.  Numerically  integrate  eq.  7,  8,  10  and  B-2  using  a  3-point  Runge-Kutta 
algorithm.  At  each  increment  solve  eq.  17  or  18  to  determine  the  optimal  re¬ 
ciprocal  temperature,  y. 

3.  At  m=m*  compare  X  to  the  desired  value,  X  *,  defining  the  error  as 

n  n 

(X  *  -  X  )2. 
n  n 

l*.  Use  a  Fibonacci  search  technique  [27]  to  adjust  the  value  of  Xj(o). 

5.  Repeat  2-1*  until  the  error  is  less  than  some  predetermined  value. 

It  was  found  convenient  to  use  a  similar  procedure  to  determine  the  iso¬ 
thermal  policy  yielding  the  same  m*  and  X  *  as  the  optimal  policy.  The  initial 
temperature  of  the  system  was  guessed  and  adjusted  to  reduce  the  error.  Also 
the  policy  consisting  of  a  step  increase  from  TA  to  T*  which  yields  the  same 
m*  and  X^*  as  the  optimal  policies  was  calculated  for  comparison.  This  policy 
would  be  the  most  extreme  temperature  rise  possible. 

The  results  of  five  such  calculations  arc  summerized  in  Table  11.  Reaction 
conditions,  conversions,  X^  and  the  dispersion  indicies  for  the  policies  that 

maximize  D  ,  minimize  D  and  reach  m*  and  X  *  isothermal ly  and  with  a  step  rise 
n’  n  n  ' 

from  Ta  toT*  are  presented.  Runs  with  the  suffix  'A'  were  calculated  without 
the  gel  effect.  The  following  conclusions  can  be  drown  from  Table  II: 

1.  The  minimum  D  is  significantly  less  than  the  isothermal  D  when  a 
strong  gel  effect  is  present  (exs.  4  and  5)  or  when  a  moderate  gel  efrect  and 
dead  end  conditions  ex!st  simultaneously  (ex.  2).  Under  these  conditions,  iso¬ 
thermal  X  .  varies  significantly  over  the  course  of  the  reaction.  Under  conven- 
tional  conditions,  minimum  ;nd  isothermal  0n 1 s  are  atout  the  same. 

2.  The  maximum  0  is  significantly  greater  than  the  isothermal  D  for 

n  n 

all  the  examples  in  Table  II-  The  gel  effect  (compare  examples  with  the  suffix 
'A'  to  those  without  it)  and  dead  end  conditions  (compare  ex.  2  to  ex.  1  and 
ex.  5  to  ex.  4)  increase  the  maximum  0^. 


I 


1  h 


3.  The  step  rise  0^  is  significantly  greater  than  the  isothermal  for 

all  the  examples  except  ex.  5-  Thp  gel  effect  and  dead  end  conditions  generally 

1  * 

decrease  the  step  rise  0  . 

n 

4.  The  maximum  0  is  significantly  greater  than  the  D  from  a  segregated 

1 

continuous  stirred  tank  reactor  (SCSTR).  Tadmor  and  Biesenberger  [2]  found 

1 

that  the  SCSTR  gives  the  broadest  MWD  of  any  isothermal  reactor  for  chain  addition 
polymerizations.  Thus,  the  maximum  effects  of  temperature  variations  op  MWD's 

are  much  greater  than  the  maximum  effects  of  residence  tinie  distributions  or 

» 

micromixing. 

1 

5.  Maximum  D^'s  resulting  from  thermal  effects  alone  are  in  the  range 

of  those  reported  [28]  for  low  density  polyethylene.  Therefore,  non- isothermal 

reactor  conditions  may  be  at  least  partially  responsible  for  these  high  disper- 

I  } 

slties. 

The  temperature  policies  and  resulting  instantaneous  and  cumulative  number 

'  1  ' 

and  weight  average  chain  lengths  are  plotted  versus  monomer  conversion  for  examples 
1,  2  and  3  on  Figures  2,  3,  and  ^ ,  respectively.  Examples  of  the  three  possible 

I 

types  of  behavior  of  isothermal  X  .  are  shown  on  these  figures. 

n,  1 

1  ■  _ 

Figure  2  shows  the  optimal  policies  for  the  case  where  the  isothermal  X 

n,  1 

decreases  then  increases  with  time.  The  temperature  policy  that  minimizes  0^ 

decreases  then  increases  to  keep  .  constant.  The  temperature  policy  that 

maximizes  Dp  is  a  step  increase  then  a  step  decrease  in  temperature  from  TA  = 

^0°C  to  T*  =  120°C  to  Tft.  This  policy  results  in  a  bimodal  distribution  with 

one  mode  formed  at  40°C,  X  .  =  850,  and  the  other  formed  at  120°C,  X  .  1  *»0 . ‘ 

n, 1  n ,  1 

The  cumulative  number  a  erage  chain  length  is  determined  primarily  by  the  higi 
temperature,  low  molecular  weight  mode,  while  the  cumulative  weight  average  chain 
length  is  determined  primarily  by  the  low  temperature,  high  molecular  weight 
mode . 

Figure  3  shows  the  optimal  policies  for  the  case  where  the  isothermal  X  , 

n,l 
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increases  with  time.  The  minimum  Dn  policy  is  one  of  increasing  temperature 

which  keeps  X  .  constant.  The  maximum  D  policy  is  a  step  decrease  in  tempera- 

n, i  n  ;  i  ! 

ture  which  again  results  in  a  bimodal  distribution. 

1  Figure  4  shows  the  optima!  policies  for  the  case  where  the  isothermal 

X  .  decreases  with  conversion.  The  minimum  D  policy  is  one  of  decreasing 
n ,  i  n  ■'  J 

temperature  which  keeps  X  .  constant.  The  maximum  D  policy  is  a.  step  increase 

n  ^  i  n 

in. temperature  which  results  in  a  bimodal  distribution. 

The  bimodal  nature  of  the  distributions  produced  by  the  maximum  0^  policy 
and  by  the  step  increase  policy  of  example  2  are  shown  in  Figure  5-  MWD's  produced 

i  ■  , 

by  the  isothermal  and  minimum  Dn  policies  are  also  plotted  on  Figure  5- 
i  The  physical  strategy  behind  the  maximum  0  policy  is  elucidated  by  the 

I  I 

numerical  examples.  is  a  measure  pf  the  standard  deviation  or,  to  use  a 

1 

physical  analogy,  the  moment  of  inertia  of  the  MWD.  The  moment  of  inertia  of 

i 

the  MWD  Will  be  maximized  if  all  the  polymer  molecules  lie  at  the  extremes  of 

t 

the  distribution  just  as  the  moment  of  inertia  of  a  body  is  maximized  if  all  its 

weight  is  concentrated  at  its  ends.  Therefore,  the  maximum  policy  seeks  to 

i  i 

form  a  bimodal  distribution  with  the  widest  pbssible  separation  between  the  modes. 

The  physical  rational  behind  the  direction  ar.d  number  of  step  changes  can  be 

seen  with  the  aid  of  Figure  6.  Inspection  of  Figure  6  shows'  that  step  changes 

in  other  than  the  directions  indicated  by  the  solid  lines  result  in  less  than 

maximum  D  'si  The  up-slopes  of  the  isothermal  X  ,  curves  for  Cases  II  and  III 
n  i  n,  i 

on  Figure  6  increase  with  increasing  gel  or  dead  end  effects.  Therefore,  the 
separation  of  thp  modes,  when  optimal  directioned  changes  are  made,  increases, 

i 

and  the  ratio  of  maximum  D  to  isothermal  D  increases.  Conversely,  non-optimal 

n  '0 

directioned  changes  increase  the  overlap  of  the  modes  and  result  in  a  decrease 

in  the  ratio  of  step  rise  to  isothermal  D  's. 

n 


t 
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EXPERIMENTAL  EXAMPLES 

Styrene  Initiated  by  AIBN  was  polymerized  using  step  change  temperature 
policies.  The  resulting  polymers  were  analyzed  by  Gel  Permeation  Chromatography 
(GPC)  to  see  if  bimodal  distributions  could  indeed  be  formed  by  temperature 
variations  alone. 

Apparatus  and  Procedure 

Styrene  (Eastman  ^Wi65  stabilized  with  tert-buty  pyrocatechol )  was  purified 
by  washing  with  2 %  KOH  followed  by  distilled  water.  The  styrene  was  then  dried 
over  anhydrous  calcium  chloride,  distilled  under  vacuum,  and  stored  on  ice  for 
not  more  than  1  hour.  Solutions  of  cold  styrene  and  AIBN  (Eastman  #61i00  2,2' 

Azobis  (2-mcthylpropioni tri le))  were  made  up  grav imet r i ca 1 1 y .  Reaction  tubes 
(Ver  Tis  drying  ampoule,  25  nil)  were  filled  wi th  about  10  ml  of  the  cold  solution. 

The  tubes  were  degassed  under  5  mm  vacuum  by  h  cycles  of  freezing  and  thawing 
and  sealed.  Tubes  were  stored  frozen  until  used. 

Prior  to  making  a  run,  the  tubes  were  allowed  to  thaw.  During  a  run, 
the  tubes  were  agitated  by  clamping  them  to  mechanical  stirrers  immersed  in 
a  constant  (-0.2°C)  temperature  bath.  To  approximate  step  temperature  changes 
2  baths  at  different  temperatures  were  used,  and  the  tubes  were  quickly  transfered 
from  one  bath  to  the  other. 

At  the  end  of  a  run  tubes  were  quenched  by  plunging  them  into  ice.  The 
tubes  were  then  opened,  the  polymer -monomer  solutions  were  weighed,  dissolved 
in  2  volumes  of  benzene  and  the  polymer  was  precipitated  in  excess  chilled  methano'  . 
After  standing  in  a  refrigerator  overnight,  the  polymer  was  separated  by  centri¬ 
fuging  at  20,000  RPM  for  1  hour.  Conversions  were  determined  gravimetrical ly. 

Tf.e  recovered  polymer  was  then  analyzed  using  a  Waters  Model  200  GPC.  Five 

Styrogel  (Waters  Associates)  packed  columns  with  nominal  pore  diameters  of  7X10^, 

k  k  3  0 

3X10  ,  10  ,  250  and  10  A  were  used.  The  ODCB  plate  count  for  the  5  column 
system  was  650  plates/ft.  The  GPC  was  calibrated  using  mono-disperse  styrene 
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EXPERIMENTAL  EXAMPLES 

Styrene  Initiated  by  AIBN  was  polymerized  using  step  change  temperature 
policies.  The  resulting  polymers  were  analyzed  by  Gel  Permeation  Chromatography 
(GPC)  to  see  if  bimodal  distributions  could  indeed  be  formed  by  temperature 
variations  alone. 

Apparatus  and  Procedure 

Styrene  (Eastman  H\kb5  stabilized  with  tert-buty  pyrocatechol )  was  purified 
by  washing  with  2%  KOH  followed  by  distilled  water.  The  styrene  was  then  dried 
over  anhydrous  calcium  chloride,  distilled  under  vacuum,  and  stored  on  ice  for 
not  more  than  1  hour.  Solutions  of  cold  styrene  and  AIBN  (Eastman  H6l\00  2,2' 

Azobis  (2-methyl propioni tri le)) were  made  up  gravimetrical ly.  Reaction  tubes 
(Ver  Tis  drying  ampoule,  25  ml)  were  filled  with  about  10  ml  of  the  cold  solution. 

The  tubes  were  degassed  under  5  mm  vacuum  by  A  cycles  of  freezing  and  thawing 
and  sealed.  Tubes  were  stored  frozen  until  used. 

Prior  to  making  a  run,  the  tubes  were  allowed  to  thaw.  During  a  run, 
the  tubes  were  agitated  by  clamping  them  to  mechanical  stirrers  immersed  in 
a  constant  (^0.2°C)  temperature  bath.  To  approximate  step  temperature  changes 
2  baths  at  different  temperatures  were  used,  and  the  tubes  were  quickly  transfered 
from  one  bath  to  the  other. 

At  the  end  of  a  run  tubes  were  quenched  by  plunging  them  into  ice.  The 
tubes  were  t.h»  n  opened,  the  polymer -monomer  solutions  were  weighed,  dissolved 
in  2  volumes  of  benzene  and  the  polymer  was  precipitated  in  excess  chilled  methano' . 
After  standing  in  a  refrigerator  overnight,  the  polymer  was  separated  by  centri¬ 
fuging  at  20,000  RPM  for  1  beer  .  Conversions  were  determined  gravimetr ical ly. 

The  recovered  polymer  was  then  analyzed  using  a  Waters  Model  200  GPC.  Five 

Styrogel  (Waters  Associates)  packed  columns  with  nominal  pore  diameters  of  7X10"’, 

4  t,  3  0 

3X10  ,  10  ,  250  and  10  A  were  used.  The  0DCB  plate  count  for  the  5  column 
system  was  650  plates/ft.  The  GPC  was  calibrated  using  mono-disperse  styrene 
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SJ- 


saniples  (Pressure  Chemicals). 

Resul ts 

Initially  four  isothermal  runs  were  made  to  test  the  validity  of  the  kinetic 

model  and  constants  (Table  l).  The  agreement  was  quite  good  so  two  sets  of 

maximum  D  and  isothermal  runs,  each  set  having  the  same  m*  and  X  *,  were  made, 
n  n 

One  set  was  made  under  conditions  requiring  a  step  decrease  in  temperature  and 

the  other  under  conditions  requiring  a  step  increase  }'  temperature.  The  results 

of  these  8  runs  are  summarized  in  Table  III. 

The  initial  isothermal  runs  are  runs  I  to  A.  Runs  5  and  6  form  one  set 

of  isothermal  and  maximum  D  runs  with  m*  =  0.25  and  X  *  =550.  The  conditions 

n  n 

were  dead  end,  so  the  maximum  policy  is  a  step  decrease  0*  =  kO°C ,  T-  =  100°C) . 

Runs  7  and  8  form  the  second  set  with  m*  =  0.33  and  Xn*  =  150.  The  conditions 

were  conventional  so  the  maximum  0  policy  ?s  a  step  increase. 

n 

For  all  8  runs  the  agreement  between  the  experimental  and  calculated  con¬ 
versions  is  excellent.  The  values  agree  within  about  10%  for  all  runs  except 
runs  2  and  6,  where  the  agreement  is  within  about  20%. 

Agreement  betv/een  the  experimental  (determined  by  GPC  and  reported  uncorrected 
for  diffusional  effects)  and  the  calculated  weight  average  chain  lengths  is  not 
as  good.  Except  for  run  6,  the  calculated  values  average  about  25%  lower  than 
the  experimental  values.  Part  of  this  is  due  to  diffusional  effects  in  the  GPC 
columns.  Experimentally  determined  average  chain  lengths  arc  compared  to  the 
values  reported  by  the  National  Bureau  of  Standards  for  their  sample  NBS-7J6 
in  Table  III.  Their  value  for  Xw  is  about  10%  lower  than  the  experimentally 
determined  value.  The  method  of  lshige,  et  al .  (29)  to  correct  chromatograms  for 
diffusional  effects  was  tried.  It  was  not  successful,  probably  due  to  skewing 
in  our  GPC.  Nevertheless,  the  agreement  between  uncorrected  experimental  and 
calculated  chain  lengths  is  sufficiently  good  to  show  that  the  maximum  0^  policy 
results  in  a  significant  increase  in  Dn  compared  to  the  isothermal  policies. 
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Normalized  chromatograms  for  runs  5  and  6,  and  7  and  8  are  compared  to 
each  other  and  to  theoretical  chromatograms  obtained  from  the  calculated  MWD’s 
for  these  runs  on  Figure  7-  Runs  5  and  7,  isothermal,  are  unlmodal,  while  runs 
6  and  8,  maximum  D^,  are  bi -modal,  as  predicted. 

To  summarize,  distributions  that  were  produced  under  non- Isothermal 
conditions  predicted  to  maximize  0^  had  dispersion  indicies  that  were  approximatly 
3  times  the  maximum  0^  that  would  be  predicted- [2]  for  an  isothermal  system. 

IDEALIZED  BATCH  REACTORS 

Having  bounded  the  problem  by  considering  temperature  variations  that  minimize 
and  maximize  the  breadth  of  the  MWD,  it  is  fair  to  ask  what  effect  actual  tempera¬ 
ture  variations  that  might  occur  ?n  polymerization  reactors  would  have  cn  the 
breadth  of  the  MWD.  Would  actual  temperature  variations  result  in  D_*s  near 
those  predicted  by  maximum  Dn  policies  or  would  they  be  closer  to  those  predicted 
for  isothermal  systems?  As  a  first  approximation  to  answering  this  question, 
temperature  variations  in  an  Idealized  batch  reactor  were  considered. 

The  idealized  reactor  is  assumed  to  be  a  batch  reactor  with  a  heat  transfer 
surface  (jacket  or  coil).  The  reactor  is  well  agitated  so  it  is  homogeneous  with 
respect  to  composition  and  temperature.  The  overall  heat  transfer  coefficient 
of  the  jacket  or  coil  is  assumed  to  be  constant,  independent  of  conversion  or 
temperature.  With  these  assumptions  the  heat  balance  on  the  reactor  can  be  written 
as: 

Ji-  „.ai  *  "ad  dm_  (22) 

<VV  dt 

where  a  is  the  reactor  diffusivity  and  N  ^  is  the  adiabatic  temperature  rise 
for  the  polymer.  Eq.  22  was  solved  numerically  in  conjunction  with  eqs.  7,  8, 
and  10  to  determine  the  effect  of  naturally  occuring  temperature  variations  on 


MWD. 


Some  temperature  profiles  for  the  bulk  polymerization  of  styrene  initiated 
by  AIBN  (lQ  =  0.3*i8  mol. /I.)  with  h(l-c)  =  1-c  and  an  initial  and  coolant  temp¬ 
erature  of  *»0°C ,  and  a  final  temperature  of  120°C  are  shown  on  Figure  8.  a 
was  varied  from  0  to  0.01  to  simulate  temperature  profiles  from  adiabatic  to 
isothermal.  Dispersion  indices  for  $  =  1  or  m  =  0.9  are  shown  on  the  curves. 

For  this  case  and  others  not  shown,  the  highest  Dn  is  r.ot  produced  by  the  adia¬ 
batic  temperature  rise  (a  =  0) .  Conditions  under  which  some  long  chains  are 
initially  produced  at  low  temperatures,  followed  by  a  rapid  temperature  rise 
producing  short  chains,  results  in  the  highest  Dn's.  , 

The  broadest  MWD  formed  under  the  conditions  of  Figure  8  (a  =  0.00025) 
is  compared  to  the  MV/D 1 s  for  the  maximum  and  minimum  policies  resulting  in 
the  same  m*,  0.261,  and  Xn*,  71,  on  Figure  The  dispersion  indices  for  the 
three  cases  are  3-**5,  7.61  and  1 . f»9 ,  respectively.  The  polymerization  conditions 
are  conventional.  The  Isothermal  case  has  approximately  the  same  MWD  and  Dn 
as  the  minimum  case.  The  temperature  variations  in  an  idealized  batch  reactor 
resulted  in  a  D  about  double  that  of  the  Isothermal  case  and  about  half  that 


of  the  maximum  0^  case.  Therefore,  conditions  exist  where  the  Dn  produced  by 
naturally  oc.currlng  temperature  variations  approaches  the  maximum  D^.  In  general 
such  conditions  will  lead  to  conventional  polymerization  where  the  isothermal 
X  ,  decreases  with  time  For  conventional  polymerizations  the  maximum  D  policy 
is  a  step  increase  in  temperature,  and  this  is  simulated,  to  some  extent,  by  the 
temperature  rise  in  the  idealized  batch  reactor.  On  the  other  hand,  the  minimum 
0^  policy  for  dead  end  polymerizations  or  polymerizations  with  a  strong  gel  effect 
is  one  of  increasing  temperature,  am  in  these  cases  the  temperature  rise  in 
the  idealized  batch  reactor  may  not  yield  a  polymer  v<ith  a  broader  MWD  than  the 
isothermal  case. 


In  *‘eal  reactors  another  factor  must  be  considered,  non-homogeneity.  Gra 
dients  and  hot  spots  In  batch  and  continuous  stirred  tank  reactois,  and  radial 
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as  well  as  axial  gradients  in  tubular  reactors  should  result  in  much  broader 
MWD's  than  those  predicted  by  the  idealized  model  used  here  [3].  Therefore, 
the  maximum  policy  may  be  a  valid  upper  limit  for  the  breadth  of  the  MWD 
from  real  reactors. 

CONCLUSIONS 

1.  The  temperature  variations  that  minimize  the  breadth  of  the  MWD  keep 
instantaneous  number  average  chain  length  constant.  Such  variations  keep  the 
MWD  distribution  from  changing  due  to  changes  in  monomer  and  initiator  concen¬ 
trations.  Such  variations  are  decreasing  temperatures  under  conventional  con¬ 
ditions  where  the  isothermal  .  decreases  with  time,  and  increasing  tempera¬ 
tures  under  strong  gel  or  dead  end  conditions  where  the  isothermal  X  .  Increases 

n>  1 

with  time. 

2.  The  temperature  variations  that  maximize  the  breadth  of  the  MWD  are 
step  changes  in  temperature  between  the  minimum  and  maximum  allowable  tempera¬ 
tures  in  the  system.  Such  variations  produce  a  bimodal  distribution  with  the 
maximum  separation  between  modes.  For  conventional  polymerization,  isothermal 

X_  .  decreases  with  time,  the  optimal  variation  is  one  step  increase  in  temperature, 
r. ,  1 

For  dead  end  polymerization,  isothermal  X  .  increases  with  time,  the  optimal 

n,  1 

variation  is  one  step  decrease  in  temperature.  For  polymerization  with  a  gel 
effect,  isothermal  X  ,  decreases  then  Increases  with  time,  the  optimal  variation 
is  a  step  increase  followed  by  a  step  decrease  in  temperature. 

3.  Temperature  variations  have  a  much  greater  effect  on  the  breadth  of  MWD's 
for  chain  addition  polymerizations  than  residence  time  distributions  or  micomixing. 

1*.  Rising  temperatures  which  may  occur  in  real  reactors  can  lead  to  a 
significant  increase  in  the  breadth  of  the  MWD  compared  to  isothermal  operation. 

In  general  this  effect  will  be  greater  for  conventional  polymerizations  than 
for  dead  end  polymerizations  or  polymerizations  with  a  strong  gel  effect. 
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NOTATION 


a 

A 


Ad,Ap,At 


c 

c* 


Ed’VEt 


f 

F 

9 

G 

AH 

H 

h 

I 


m 

m* 

M 

M 

) 


initial  reactant  ration,  fln/M 

o  o 

heat  transfer  area 

frequency  factor  for  initiator  decomposition,  propogation,  and 
termination,  respectively. 

initiator  conversion,  (1  -  i / I Q) 

target  initiator  conversion 

heat  capacity 

dispersion  index,  X  /X 
r  w  n 

activation  energies  .'or  initiator  decompos i t ion, propogat ion? 
and  termination. 

initiator  efficiency 

A-7 

eq.  A-7 

(5-c)n 

initiator  concentration 
initial  initiator  concentration 

kinetic  constants  for  initiator  decomposition  and  for  propogation 
respectively 

2 

termination  constant,  (k  +  k^)  *=  k^fyj-g  (m) 

kP(2f,okd/kto)i  =  Arixp('Eiy) 

'on',kd  =  A2  exp  ('e2y) 

M  k  2/k  =  A.  exp(-E.y) 

o  p  t  5  i 

monomer  conversion,  (1-M/Mo) 
target  monomer  conversion 
monomer  concentration 


moment  generating  function  defined  by 

gel  effect  function,  g(m)  = 

moment  generating  function  defined  by 

heat  of  polymerization 

u  ...  .  I  X.  tdX./dt) 

Hamiltonian,  ,  I  i 


Initial  monomer  concentration 
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ad 


R 

R 

9 

t 

T 

Tf 

T 

o 

V 

U 

x 

x 


x  * 

n 

X  . 
n,i 


y 

Y* 


Greek  letters 


a 

T) 


concentration  of  dead  polymer  of  length  x 

adiabatic  temperature  rise,  -AH/(MW*cv) 

concentration  ol*  growing  chains  of  length  x 

n^  moment  of  the  dead  polymer  distribution,  defined  by  eq. 

A-8.  ' 

concentration  of  primary  radicals 

gas  constant 

time 

absolute  temperature 
final  temperature 

I 

initial  temperature,  coolant  temperature 
Reactor  volume  , 

Overall  heat  transfer  coefficient 

state  vector  with  components  Xj 

number  average  chain  length 

target  number  average  chain  length 

Instantaneous  number  average  chain  length,  dm/dCo 

weight  average  chain  length 

reciprocal  temperature,  l/R^T 

maximum  y,  corresponds  to  physically  constrained  minimum  temperature. 

V 

minimum  y,  corresponds  to  physically  constrained  maximum  tempera¬ 
ture,  T* 

stationary  temperature  policy.  Solution  to  3H/3y  =  0 

moment  of  the  living  polymer  distribution,  defined  by  eq.  A-8 

reactor  d i f fus i v i ty , (UA/Vp  cv) 

ratio  of  k^/k^ 

adjoint  variable,  -3H/3X. 
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dimensionless  temperature,  (T-T^J/CT^-T  ) 
density 

standard  deviatioi. 

dimensionless  moment  of  the  dead  polymer  distribution,  q  /M 

no 
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APPENDIX  A.  DERIVATION  OF  STATE  EQUATIONS 

With  the  assumption  of  constant  density  and  that  primary  radicals  do 
not  terminate  differential  mass  balances  for  a  batch  reactor  can  be  written 
from  reactions  1-5. 


dj_  _ 
dt 

-kdh  (1) 

dM 

03 

-k  M  E 

p  - 

k.HR 

dt 

P  n=l 

n 

i 

dR 

dt 

-  2fkdh(D 

-k.HR 

i 

00 

dt 

=  k  .MR  - 

i 

k  MP.  - 
P  1 

(ktc  + 

ktd)  P,  S  .  P 

1  n=i  n 

dP 

00 

X 

dt 

=  k  M  (P 
p  x- 

1  -  Px> 

'  (ktc 

*  k«’  p>  A  p" 

dM 

_ X 

X 

=  ik 

-i 

E  P.l 

3  .  + 

00 

k.jP  E  P. 

dt 

tc  . 

=1  1 

x-  1 

td  x  n=)  n 

(A-l) 
(A-2) 
(A- 3) 


(A-l.) 


(A-5) 


(A-6) 


The  moments  of 

130] . 

00 

G(r,x)  =  E 
x=l 


the  MWD  are  derived  by  means 


CO 


rXPx  i  F(r,x)  =  Z 

x=l 


of  moment  generating  functions 
r\  (A-7) 


from  which  the  moments  of  the  living  and  dead  chain  distributions  can  be 
obtained  by  differentiation 


9  nG  ( I ,x)  ;  qn  =  9  PF  ( 1 ,x) 

arn  8  rn 


(A-8) 


X  X 

G  and  F  are  derived  by  summing  r  times  eq.  A-5  and  r  times  eq.  A-6,  respectively. 


dG 

dt 


k  .HR 


k  H  E  rXP  +  k  M  E  rX+,P 


x=l 


x=  I 
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the'QSSA  to  A*3,  -11  to  -13  yields: 


2 


O 


/  ?fkdh{l)\i 

\ktc+ktd  / 


k  M 


k  +k^ 
tc  td 


22  '  2\HTI 


(A 


Substitution  of  A-18  and  the  deffinitions  of  m,  c,  n,  the  dimensionless 

,  i 

moments,  C j »  and  the  lumped  kinetic  constants,  k j ,  k 2  and  k^,  in  A- 1 ^  to 
-16  yields  the  state  eqi.  1 0- 13- 


APPEND I X  B.  DERIVATION  OF  CRITERIA  FOR  THE  MINIMUM  D  POLICY, 
- n* — - - 

,  When  n  is  constant  the  state  of  the  system  is  described  by  equations  7,8 

I  , 

and  10.  The  optimal  policy  makes  the  Hamiltonian  a  maximum  and  zero.  The 
Hamiltonian  is  given  by: 

H  •=  0  =  A,  [kjh^O-c).  (l"m)  +  2ak2h(l-c)]  +  h(!l-c) 

9  2 

+  A,  f  (2+t,)k_  -^4-  +  2(l+n)k.hi(l-c) 

i  3  92  1  9 

+  2nak2h(l-c) ]  (8-1) 


When  f  is  constant,  the  adjoint  variables  (Aj  =  -3H/3X 


dA, 


dt 


[A,  +  2A3(l+n^]  k,h  dm  t  g 


id  f  Ijnn 

’g 
2 


.  ,  _  , ,  d  (l-m)‘ 

3(2  n)k3  dm  02 


dX2 

dT 


[A,  +  2A3(l+n)]  k, 


(1-m)  d 
g  dc 


(h‘) 


•  [ 2a A  j  +  A2  +  2anA3]  k2 


(h) 


I )  are  described  by: 


(8-2) 


(6-3) 

(B-k) 
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The  boundary  conditions  on  the  state  eqs  7,  8  and  10  and  adjoint  eqs.,  B-l 
to  -A,  for  the  case  where  polymer  is  not  initially  present  are: 


m  = 


c  = 


X. 


t  =  0 
0 

0 


t  =  G 


c*  = 


min 


(2-n)  m* 

aX  * 
n 


(B-5) 


Therefore,  X^  =  -l 

Along  the  optimal  path,  where  H  =  0,  the  3H/3y  is  given  by: 


=  (E2-E,)  {x,-2-2n)klhi-iija)-  ♦  (E3-E2)  (2*n)k3  -(^ 


(B-6) 


From  the  definition  of  kj ,  k^  and  k^  the  following  identities  can  be  proven: 


2ak„ 


(E2-E,)  »  -2(E3-E2) 


Substitution  of  (B-7)  into  (B-6)  yields: 
—  =  (E2-E,)  (X | -2-2n) k| h  — - 


2 (e  -e  )  (2tPl 
/U2  r  2a 


(B-7) 


kl  (1-m)2 


(B-8) 


Along  the  stationary  policy,  y*  <  y  $  y*,  DH/3y  = 
32H  ,,  ,2  (2+n)  kl  (1-m)2 

^  -  -(e2-Ej)  — 2T"  k2  — r~ 


0: 

<  0 


(B-9) 


Therefore,  the  stationary  policy  can  form  part  of  the  optimal  path. 

Applying  the  criteria  for  the  optimal  path,  eq.  16,  to  and  rearranging 


B-8  yields: 
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k,  ( 1  —ns) 
akjh^g 


i  y  *  y* 

-  (x.  -2-2n)  ;  y  =  ys 

(2+n)  '  y  =  r* 


(B-IO) 


Substitution  of  eq.  13  into  B-10  yields: 

<  (x\ ' 2  ~ 2r))  2 y =  y* 

n,i  >  (2  +  n)  (2  -  n)  T2~-  nT  '  J  =  \ 


(B* 1 1 ) 


But  along  the  stationaiy  path,  3H/3y  =  0.  Substituting  eq.  B-10  into  B-2  yields 

f,  i  -  ,  i  ^  ,l-m  . 

dt  *'  =  -t*j~2-2n]  kjh  (  g  ) 


■<E2-E1> 

2(VE|> 


(X.-2-2.,)  k,  - [  2  J  -4-  (-J3— ) 

1  1  (l-m)  9  dm  9 


a  o  ;  X,  =  const. 

Q  l  I 


(B-12) 


And  therefore  along  y 


X  .  =-  const, 

n,  i 


(B-13) 


as  predicted  by  Hoffman,  et  al.  [6]. 

APPENDIX  C.  DERIVATION  OF  CRITERIA  FOR  MAXIMUM  D  POLICY 

-  - - - _n - 

The  state  of  the  system  is  described  by  eqs.  7>  8  and  10.  The  Hamiltonian 
is  given  by  eq.  B-l,  and  the  adjoint  variables  are  described  by  eqs.  B-2  to  B-4. 
The  boundary  conditions  for  the  state  and  adjoint  variables  with  the  exception 
of  X^  are  given  by  B-5.  The  boundary  conditions  for  X^  are: 


t  *=  0 


t  =  0 


Therefore,  X^  =  +1 . 

Along  the  optimal  path  where  H  =  0: 


|  -  (Ej-E,)  (X,  *  2<2,)k,hi  iii>  .  2(E2-E|) 
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Along  the  stationary  policy,  y*  <  y  4  y.,.,  3H/3y  =  0: 


<VEi>' 


Therefore,  the  stationary  policy  minimizes  rather  than  maximizes  the  Hamiltonian, 
and,  from  the  Maximum  Principle,  cannot  form  part  of  the  optimal  path. 

Rearrangement  of  C-2  and  the  application  of  eqs.  3»  6  and  C-3  yields  the 


criteria  for  the  optimal  policy: 
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